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ABSTRACT 

We describe the measurement of the beam profiles and window functions for the Atacama Cosmology Telescope 
(ACT), which operated from 2007 to 2010 with kilopixel bolometer arrays centered at 148, 218, and 277 GHz. Maps 
of Saturn are used to measure the beam shape in each array and for each season of observations. Radial profiles 
are transformed to Fourier space in a way that preserves the spatial correlations in the beam uncertainty to derive 
window functions relevant for angular power spectrum analysis. Several corrections are applied to the resulting 
beam transforms, including an empirical correction measured from the final cosmic microwave background (CMB) 
survey maps to account for the effects of mild pointing variation and alignment errors. Observations of Uranus 
made regularly throughout each observing season are used to measure the effects of atmospheric opacity and to 
monitor deviations in telescope focus over the season. Using the WMAP - based calibration of the ACT maps to 
the CMB blackbody, we obtain precise measurements of the brightness temperatures of the Uranus and Saturn 
disks at effective frequencies of 149 and 219 GHz. For Uranus we obtain thermodynamic brightness temperatures 
7^ 49 = 106.7 zb 2.2 K and 7^ 19 = 100. 1 ± 3. 1 K. For Saturn, we model the effects of the ring opacity and emission 
using a simple model and obtain resulting (unobscured) disk temperatures of T s 149 = 137.3 =b 3.2 K and T s 219 = 
137.3 ±4.7 K. 

Key words : cosmology: observations - planets and satellites: individual (Saturn, Uranus) 

Online-only material: color figures 


1. INTRODUCTION 

The Atacama Cosmology Telescope (ACT) is a 6 m, off- 
axis telescope designed for millimeter wavelength observations 
of the cosmic microwave background (CMB) with arcminute 
resolution (Fowler et al. 2007). The telescope is located in 
Northern Chile, at an elevation of 5190 m, where atmospheric 
conditions are excellent for microwave observations from April 
through December. From 2008 to 2010, ACT observed the sky 
with bolometer arrays operating in frequency bands centered at 
148 GHz, 218 GHz and 277 GHz. The choice of bands provides 
the ability, in principle, to spectrally discriminate between signal 
from the CMB, the Sunyaev-Zel’dovich effect (Sunyaev & 
Zel’dovich 1970), and dusty sources. The 148 GHz array was 
also in operation during a short observing season at the end of 
2007. 


The ACT scientific program has produced measurements 
of the CMB angular power spectrum at multipoles spanning 
i ~ 300-10,000 (Das et al. 2013; Dunkley et al. 2013), probing 
the primary and secondary anisotropies (Sievers et al. 2013). 
ACT has produced source catalogs (Marriage et al. 201 lb), and 
galaxy cluster samples reaching to high redshifts (Marriage et al. 
2011a; Hasselfield et al. 2013). 

The telescope beam acts as a low-pass spatial filter on the 
astrophysical signal from the sky. The effects of this filtering 
are relevant, to varying degrees, to all ACT science results. 
Interpretation of the angular power spectrum, in particular, relies 
on its decomposition into components that are correlated across 
large ranges in angular scale. The extraction of cosmological 
parameters associated with such components depends on an 
understanding of the beam shape and the covariant features of 
the window function uncertainty. 
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In this work we describe the processing of planet observations 
to obtain the telescope beams, and thus the instrument window 
function, for each detector array for each season of operation. 
We also discuss the use of Uranus observations to monitor beam 
shape and calibration. Because the ACT maps are ultimately 
calibrated to the CMB blackbody via the dipole calibration 
of Wilkinson Microwave Anisotropy Probe ( WMAP ; Jarosik 
et al. 2011), we obtain precise measurements of the Uranus 
and Saturn disk temperatures without resorting to interpolation 
across unknown source spectra or other calibration transfer 
standards. 

A detailed description of the instrument can be found in S wetz 
et al. (2011), and the main data reduction and map-making 
pipeline is described in Diinner et al. (2013). Hincks et al. (2010) 
present an analysis of the ACT beams from the 2008 season; the 
present approach extends this earlier analysis. 

In Section 2 of this paper we discuss the planet observations 
and maps. In Section 3 we obtain the telescope window function 
through an analysis of the planet maps and the application of 
corrections for several systematic effects. In Section 4, we 
interpret the apparent brightness of Uranus to calibrate the 
sensitivity of the telescope to atmospheric water vapor and focus 
variation. In Section 5 we use the CMB -based calibration of 
the ACT maps to obtain absolutely calibrated measurements of 
Uranus and Saturn brightness temperatures. 

2. OBSERVATIONS AND MAPPING 
2.7. Observations 

The ACT observation strategy is to scan over small ranges of 
azimuth angle while keeping the boresight altitude fixed at ~50° 
to minimize systematic effects due to altitude variation. The 
CMB fields are observed as they rise and set through two central 
azimuth pointings. Planet observations are made at this same 
telescope altitude angle, by briefly interrupting CMB survey 
scans and re-pointing in azimuth to the planet location. 

ACT observed a planet approximately once per night of 
operation in 2007-2010, but only high quality observations 
of Saturn and Uranus are considered in this analysis. Saturn 
and Uranus were preferentially targeted, as Saturn’s brightness 
makes it useful for beam profile measurements, and Uranus’ 
stable brightness is a convenient calibration source. While 
Uranus was visible throughout each observing season, Saturn 
was available only at the beginning and end of each season. 
Jupiter and Mars were, respectively, too bright and too low in 
the sky to be useful for calibration work. In Table 1 we list the 
number of high quality Uranus and Saturn observations achieved 
in each season. 

2.2. Planet Maps 

As described in Diinner et al. (2013), maps are made from the 
ACT time stream data according to a maximum likelihood tech- 
nique in which the noise is described in the time domain while 
the signal is assumed to be spatially coherent (Tegmark 1997). 
The bolometer time stream data for planet observations are 
pre-processed in a similar fashion to the CMB survey data. The 
time stream data are calibrated, based on detector load curves, 
after deconvolving the effects of the detector time constants and 
low-pass readout filtering. Detectors are then screened for qual- 
ity based on their projections onto a common-mode computed 
using a fiducial set of well-behaved detectors. Detector screen- 
ing is performed with the planetary signal masked out. Maps 


Table 1 

Number of Saturn and Uranus Observations Selected for 
Calibration Purposes for Each Season and Array 



2007 

2008 

2009 

2010 

Saturn 

148 GHz 

22 

22 

13 

35 

218 GHz 


23 

9 

28 

277 GHz 


16 

5 

23 

Uranus 

148 GHz 

16 

37 

94 

113 

218 GHz 


33 

75 

102 

277 GHz 


21 

16 

74 


Notes. While Uranus was visible throughout each observing 
season, Saturn was available only near the beginning and end 
of each season. The number of successful planet observations is 
lower at higher frequencies because of increased sensitivity in to 
atmospheric contamination these arrays. 

are made using a dedicated code that is optimized for high res- 
olution mapping of single observations; we have confirmed that 
the main season map pipeline code produces compatible results 
when run on the same data. Each map is solved iteratively, con- 
verging after fewer than 10 iterations. On long time scales the 
bolometer data are dominated by common-mode signal from 
the atmosphere. The map-maker includes a noise model that 
strongly suppresses common-mode components at frequencies 
below 0.2 Hz, corresponding to angular scales of approximately 
3°, which is larger than the array size. We assess the impact of 
this mild high-pass filtering on our beam measurement at large 
angular scales in Section 3.4.1. 

2.3. Map Selection 

While planets were observed almost every night, observations 
do not always yield successful maps. The primary reason 
for failed mapping is a low number of live detectors during 
periods of high sky loading. This is a greater problem at higher 
frequencies, and for the 277 GHz array in particular leads to 
many maps with incomplete sampling of the planet signal; such 
maps are discarded. 

Planet observations are sometimes made after sunrise. While 
such observations allow us to study the magnitude of pointing 
and beam focus changes due to thermal deformation of the 
telescope structure, they are excluded from the present analysis. 
The sunrise cut is also applied to CMB observations. 

A small number of maps are cut due to having substantially 
higher noise levels than other maps of the same planet for the 
same array and season. The remaining maps are used to establish 
the telescope beam and calibration parameters, as described in 
the following sections. Season mean beam maps, obtained by 
aligning and averaging the selected Saturn maps for each season 
and array, are shown in Figure 1 . The Airy pattern is easily seen 
in the maps at 148 GHz and 218 GHz. The horizontal streaks 
are parallel to the telescope scan direction. 

3. BEAMS 

In this section we obtain measurements of the telescope beams 
for each season and array. Saturn is bright and much smaller 
than the beam solid angle, permitting the characterization of 
our beams with high signal-to-noise ratio. Saturn maps are re- 
duced to radial profiles, which are then modeled with appropri- 
ately chosen basis functions and transformed to Fourier space. 
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Figure 1. Beam maps for each season and array, formed from the average of selected Saturn maps. Gray scale is logarithmic. Each panel is 12' by 12'. Each contributing 
Saturn map incorporates all live detectors in the array. The well-resolved Airy rings indicates that the relative detector offsets are accurately known. The frequency 
labels correspond to the effective frequencies for CMB spectrum radiation (Swetz et al. 2011). 


The resulting transform, the beam modulation transfer function, 
is corrected for a number of systematic effects to produce a 
window function and i - space covariance suitable for use with 
the ACT CMB survey maps. 

3.1. Radial Profiles 

While the instantaneous ACT beams are slightly (<10%) 
elliptical in cross-section, we analyze them as though they 
were circular, and work exclusively with the “symmetrized” 
(i.e., azimuthally averaged) radial profiles. The resulting beam 
measurements are suitable for any analysis that incorporates 
a similar simple azimuthal average of the signal under study. 
In particular, the cross-linking of CMB survey observations 
is such that the telescope beam contributes to the map, with 
roughly equal weight, at two orientations that differ in rotation 
by approximately 90°; the net effective beam is very nearly 
circular. The goal of our analysis is thus to characterize the 
symmetrized beam. 

Radial profiles are obtained for each Saturn observation by 
averaging map data in radial annuli of width 9", producing 
average profile measurements y t at radii Of . On scales smaller 
than a few arcminutes, the covariant noise is sub-dominant to 
the planetary signal and thus the shape of the core beam is very 
well measured by even a single observation of Saturn. However, 
at larger angular scales there is non-negligible contamination 
from the atmosphere. This means that the asymptotic behavior 
of the beam cannot be separated from the fluctuating background 
without making some assumptions about the beam behavior far 
from the main lobe. 


The illumination of the ACT optics is controlled by a 
cryogenic Lyot stop at an image of the circular exit pupil; 
this leads to a beam shape described by an Airy pattern (for 
monochromatic radiation) and to the expectation that the beam 
pattern (including the effects of finite detector size and spectral 
response) will decay asymptotically as I/O 3 , where 0 is the 
angle from the beam peak. We fit this model to the radial profile 
data in order to fix the background level of the map and to 
provide a model for extrapolating the beam to large angles. 

The background level of the map and the large-angle behavior 
of the beam are measured by fitting the model y{6) = A x 
(l'/0) 3 + c to the binned profile data in the range 6a < ft < 0b, 
with the range chosen such that the profile measurements have 
fallen to below 1% of the peak but are still signal dominated. 
The points yt with ft < 6a, corrected for the background level c 
and renormalized so that y(0 = 0) ~ 1, constitute the “core” of 
the beam profile. The radial profile beyond 6a is referred to as 
the “wing.” The fit ranges used for each array, and the amplitude 
of the beam at 6a relative to the peak, are given in Table 2. An 
example of the wing fit for each array is shown in Figure 2. 

For each array, the multiple observations of Saturn in each 
season are combined by taking the mean of their core profile 
points and wings. The season mean beam profile is thus 
described by points y* (for each ft < 6a) and a mean wing 
fit parameter A. The full covariance matrix of {y \ , . . . , y n , A} is 
computed from the ensemble of Saturn profiles and, as described 
in the next section, is used to propagate covariant error on large 
angular scales into the i - space beam covariance matrix. 

The solid angle and FWHM of the beam (as well as the 
approximate ellipticity of the non- symmetrized telescope beam) 
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Table 2 

Summary of Beam Parameters by Array and Season 


Array 

Season 


Wing Fit (Section 3.1) 

Transform Fit (Section 3.2) 

Beam Properties 


Hincks et al. 

Solid Angle 
(1CT 9 sr) 

6a 

6b 

B(9 a ) 

(io- 3 ) 

f max 

^rnodc 

Solid Angle 
(1(T 9 sr) 

FWHM 

(arcmin) 

Ellip. 

148 GHz 

2007 

4' 

8' 

2.40 ± 0.05 

21900 

11 

224.5 ± 1.4 

1.364' 

1.098 



2008 

4' 

8' 

2.03 ± 0.04 

18400 

9 

216.6 ± 1.6 

1.373' 

1.065 

218.2 ±4.0 


2009 

4' 

8' 

2.15 ±0.09 

16400 

8 

215.8 ±2.6 

1.378' 

1.039 



2010 

4' 

8' 

2.13 ±0.04 

15800 

9 

215.4 ± 1.8 

1.381' 

1.033 


218 GHz 

2008 

2' 

4' 

6.54 ±0.12 

25500 

6 

116.5 ± 1.1 

1.092' 

1.021 

118.2 ±3.0 


2009 

2! 

4' 

7.07 ± 0.29 

25800 

6 

123.9 ± 1.7 

1.057' 

1.026 



2010 

2! 

4 ' 

7.34 ± 0.21 

25100 

6 

125.0 ± 1.8 

1.015' 

1.027 


277 GHz 

2008 

2! 

4' 

5.49 ±0.12 

28400 

7 

98.1 ± 1.9 

0.869' 

1.049 

104.2 ± 6.0 


2009 

2! 

4 ' 

5.22 ± 0.26 

32100 

7 

95.2 ± 4.0 

0.879' 

1.118 



2010 

2! 

4' 

5.65 ±0.18 

28600 

7 

98.9 ± 2.0 

0.870' 

1.103 



Notes. Radial profiles are obtained within a distance 6 a from the peak; wing and baseline parameters are fit to data over the range 6a to 6b- The level B (6a) of 
the beam, relative to the peak, at the edge of the wing fit region is shown for reference. Parameters i max and n mo de describe the basis functions used for fitting 
a beam model to the radial profile data in the beam core. The solid angle and FWHM of the symmetrized instantaneous telescope beam are provided, along 
with the ellipticity (defined as the ratio of major to minor axes of the unsymmetrized mean beam). The solid angles for the 2008 season may be compared to 
the values from the independent analysis of Hincks et al. (2010). 



Figure 2. Wing model fitted to binned radial profile data. Points shown, and the 
best-fit model (solid line) are from a typical individual Saturn observation for 
each array. A baseline is fit simultaneously and has been removed from the data 
points. The shaded area represents the standard deviation of all individual best 
fits to wing models for the 2010 season. 

are shown in Table 2. Solid angles may be compared to the values 
presented in Hincks et al. (2010). The reduction in solid angle 
uncertainty is primarily due to our choice to fit the wing over a 
range of radii closer to the beam center. We have doubled the 


error from its formal value to account for systematic variation 
in solid angle results as different fitting ranges are used. 

3.2. Harmonic Transforms and Window Functions 

Sky temperature data may be compared to predictions based 
on cosmological models through statistical measures such as 
the angular power spectrum of temperature fluctuations. The 
effect of the observing strategy and telescope beam on the 
measurement of correlation functions in CMB data has been 
discussed in, e.g., White & Srednicki (1995). For ACT, it is 
important to understand the beam and its covariance over a 
broad range of angular scales. 

A Gaussian random temperature field T(n) at position n 
on the sphere may be characterized by coefficients Q of the 
decomposition of its auto-correlation function: 

00 21 + 1 

<r(n)r(n')) = V — CtP t { n n'), (D 

tr 4jr 

where Pi are Legendre polynomials. Measurements T(n) of 
the temperature field with a given telescope and observing 
strategy will differ from the true sky temperature; for ACT the 
primary effects are an uncertain calibration (which we neglect 
for now) and the telescope beam. For the case of an azimuthally 
symmetric beam that does not vary in time or with position on the 
sky, the auto-correlation function of the observed temperature 
map T (n) will be related to the Cp describing the T (n) by (e.g., 
White & Srednicki 1995) 

<f (n)f (n')> = V n n)Bj, (2) 

i 4n 

where Bi are the multipole coefficients of the beam, defined by 
00 1(21 + 1) 

B(0) = T bJ— PA cos 6). (3) 

tr v 4?r 

In this context, Bj is referred to as the window function. 
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The ACT beams are sufficiently concentrated that we may 
work in a flat sky approximation and take 6 to be a ra- 
dial coordinate in the plane. The radial component of the 
Fourier conjugate variable will be £ . The azimuthally symmetric 
two-dimensional Fourier transform B{£) of B(6) may be used 
to obtain the multipole coefficients according to Bi = B(£) 
with entirely negligible error. We will thus work in the flat sky 
approximation, but write Bg (with i in subscript) to emphasize 
the correspondence to the multipole coefficients. 

We consider contributions from the wing and from the core 
separately. The transform of the wing (which is truncated below 
0 A ) is easily obtained. For the beam core, computation of the 
transform from the binned radial profile requires us to interpolate 
the profile data y t . We achieve this by expressing the real space 
beam as a sum of basis functions which we know to have spatial 
frequency cut-offs near those determined by the telescope optics 
in each frequency band. 

The primary mirror has a diameter D = 6 m; this size limits the 
optical response to spatial frequencies below £ max ~ 2n D/X. 
A natural basis for the Fourier space beam is then provided 
by the Zernike polynomials R^ip), where p = £/£ max . (The 
Zernike polynomials R™(p, 0) are orthogonal and complete on 
the unit disk; here we impose azimuthal symmetry and need 
only consider even n and m = 0.) We thus model the real space 
beam with basis functions f n {6 ), for non-negative integer n , 
that are proportional to inverse Fourier transforms of R^ip): 

f n (Q) — ^ 2w+1 ^^ max ) ( 4 ^ 

^^max 

In practice, for some value of £ max and a mode count n mo d e 
the beam model is expressed as 

^mocle 1 

B(0)= a n f n {0), (5) 

n— 0 

where the coefficients a n are fit parameters. The model is fitted 
to the data by minimizing the x 2 of the residuals y t — B(6t), 
accounting for the full covariance of they/ measurements. Values 
of £ max and n m 0 d e are adjusted to obtain a fit with the minimum 
number of modes necessary that gives a x 2 P er degree of 
freedom equal to unity. The fit yields coefficients a n and the 
covariance matrix of the errors, (8a n 8a n ). 

Details of the fitting parameters, including £ max , n mo d e and 
reduced x 2 for each season and array, are presented in Table 2. 
The choice of basis functions results in a satisfactory expansion 
of the beam using approximately half as many modes as data 
points. 

After fitting the coefficients a n in real space, we obtain the 
corresponding beam transform, 

B t = j2 anF ? +w t' (6) 

n 

where F” are the Fourier transforms of f n (6), truncated above 
6 At and wt is the Fourier transform of the extrapolated wing, 
from 6 a to infinity. 

The resulting beam transforms are shown in Figure 3. While 
the beams for the 2008 season are shown, the beam features are 
not substantially different between seasons. The figure shows 
the separate contributions of the core and wing to the beam 
transform. The wing contributes significantly only at low £ 
(i.e., below 1000 at 148 GHz, and below 2000 at 218 GHz 
and 277 GHz). 


While this is very similar to the procedure described in Hincks 
et al. (2010), in this work we include covariance between radial 
profile data points, and between the radial profile and the wing 
fit parameter. This results in a natural propagation of beam errors 
on all spatial scales into the fitted beam transform and its £- space 
covariance matrix. Our treatment of beam errors is discussed in 
the next section. 

3.3. Beam and Calibration Covariance 

Using the covariance matrix of the fit coefficients obtained in 
Section 3.1, we also obtain the covariance in £ space: 

= (8Bi8Bir) 

= F?Fp8a n 8a n ') + ( Sw t 8w e ,) 

n,n' 

+ Y^ F "(Sa n Swr)+y2( Sw eS an )Fl,. (7) 

n n 

This covariance includes contributions from the wing fit error 
and the correlation between the wing fit error and the coefficients 

a n . 

This covariance matrix is an £-space representation of the 
covariant features in the radial profiles that contribute to the 
Bi . However, because the calibration of the CMB survey maps 
is established at a particular angular scale, we must recast the 
beam covariance into an appropriate form. This is effectively 
the same procedure applied in Page et al. (2003). 

The absolute calibration of the ACT maps at 148 GHz is 
obtained by cross-correlation of the 2008 Southern and 2010 
Equatorial maps with the WMAP 7 year maps (Jarosik et al. 
2011) at 94 GHz of the same sky region, over angular scales 
with 300 < £ < 1 100 (Hajian et al. 2011; Das et al. 2013). This 
leads to an absolute calibration of the ACT maps centered near 
£ = 700. The 2009 season Equatorial maps at 148 GHz are then 
calibrated, over 500 < £ < 2500, to the 2010 season Equatorial 
maps. The subsequent calibration of the 218 GHz maps to the 
148 GHz maps is performed in the signal-dominated regime 
with 1000 < £ < 3000. The 277 GHz maps have not yet been 
calibrated to WMAP. For each season and array, we factor out 
the beam amplitude at effective calibration scale £ = L, where 
L = 700 for the 148 GHz array, and L = 1500 for the 218 GHz 
and 277 GHz arrays. 

For a celestial temperature signal described by multipole 
moments , we measure a map M^ m in detector power units, 
which is related to T by 

M^ m = GiTi tm , ( 8 ) 

where Gi = G$Bi is the product of the beam (normalized such 

that Bi = o = 1) and a global calibration factor Go that converts 

CMB temperature to map units (i.e., to detector power units). 

Our calibration to WMAP is a comparison of M and T at 
£ = L and is thus a measurement of G L that is independent of 
our beam uncertainty. This leads us to recast Equation (8) as 

= G L biT^ m , (9) 

where 

Bi 

= (10) 

Bl 

Since the errors in bi and G L are not correlated, the uncertainty 
in Gi is described by a covariance matrix Ta> that is the sum of 
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148 GHz 218 GHz 277 GHz 








Figure 3. Upper panels: beam transforms Bi for each array. Beams shown are for the 2008 season. The total beam is the sum of the contribution from the core of 
the beam and the extrapolated wing. Lower panels: diagonal error from the covariance matrix for the renormalized beam. The beam normalization has been fixed at 
calibration t = 700 (1500) for the 148 GHz (218 GHz and 277 GHz) array(s), as described in Section 3.3. Note that these curves show uncertainty, not systematic 
trends, and they do not include additional uncertainty from the empirical corrections of Section 3.5. The full covariance matrix shows an anti-correlation between the 
beam error at angular scales above the calibration i and beam error below the calibration i. The window function is Bf and thus its fractional error is 2 SBi/Bi. 


a calibration error from the measurement of Gl, and a term due 
to correlated error in the bf. 


2 

Tu> = bgb^{8Gi) 2 + — b{L L p — limbi' + bgbplj LL ] . 

b l 


( 11 ) 


The term in square brackets is the normalized beam covariance 
that accompanies ACT data releases. The diagonal beam error 
(i.e., the square root of the diagonal entries of the normalized 
beam covariance) is shown for each season and array in Figure 3. 
At high l, the error in the effective season beam is dominated by 
an empirical map-based correction, which is not included here 
(see Section 3.5). 

While the fitted beam and covariance are an accurate descrip- 
tion of the binned radial beam profile, they must be corrected for 
a number of systematic effects prior to being used to interpret 
ACT maps. 


3.4. Correction for Systematic s 

The beams computed in the preceding section are corrected 
for a number of systematic effects that would otherwise bias the 
resulting transforms relative to the true telescope beam. These 
are briefly described below. 

3.4.1. Mapping Transfer Function 

Because our map-making procedure includes time domain 
high-pass filtering, we might expect poor reproduction of large 


spatial scales. We study the mapping transfer function by 
injecting a simulated signal into telescope time stream data and 
comparing the output map to the input map in Fourier space. 
The transfer function deviates significantly from unity only on 
angular scales larger than 20'. The wing fit is performed at 
somewhat smaller radii than this, and we have confirmed that 
the wing fit and extrapolation reproduces the input signal, on all 
scales, to 0.1%. 

A very small (~0.1%) correction due to the 3'.'5 pixelization 
of the planet maps is applied by dividing the beam transform by 
the azimuthal average of the analytic pixel window function. 

3.4.2. Radial Binning of Planet Maps 

The binning of the planet map pixel data into annuli has a 
slight impact on the inferred beam transform. This is quantified 
by evaluating the harmonic transform of data points taken 
from a model of the radial beam profile, and comparing it 
to the harmonic transform of points that include simulated 
binning of the radial beam profile. This resulting correction 
is approximately 1% at i = 10, 000. 

3.4.3. Saturn Disk and Ring Shape 

The Saturn angular diameter of approximately 18" is large 
enough that we do not treat it as a point source. For each season 
and array, we deconvolve Saturn’s shape assuming that it is 
a disk with solid angle equal to the mean solid angle for all 
Saturn observations that contribute to the mean instantaneous 
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Table 3 

Pointing Variance and Beam Correction Parameters 


Season 

Array 

2 

°planet 

(arcsec 2 ) 

V 

(arcsec 2 ) 

Attenuation 
at l = 5000 

2007 

148 GHz 

53 ± 12 



2008 

148 GHz 

24 ±4 

25 ± 10 

0.993 ± 0.003 


218 GHz 

23 ±5 

143 ±51 

0.962 ± 0.013 


277 GHz 

11 ±2 



2009 

148 GHz 

58 ±9 

26 ±6 

0.993 ± 0.002 


218 GHz 

43 ±8 

46 ±27 

0.988 ± 0.007 


277 GHz 

49 ± 18 



2010 

148 GHz 

31 ±3 

29 ±6 

0.992 ± 0.002 


218 GHz 

18 ±2 

-5 ±23 

1.002 ± 0.006 


277 GHz 

25 ±3 




Notes. Pointing variance cr 2 lanet is measured from planet observations, 
and the season effective beam correction parameter V is measured from 
point sources in the full-season CMB survey maps, where available. The 
parameter V is used to correct the high-f beam for use with the survey 
maps; the resulting attenuation or inflation of the beam at t — 5000 
is provided for reference. The covariant error due to this correction is 
included in the season effective beam covariance matrix. 

beam measurement. This correction is approximately 2% at 

1 = 10 , 000 . 

While Saturn’s rings complicate the spatial distribution of the 
planetary signal, these effects are negligible at t < 20,000, aside 
from an overall change in average brightness, and need not be 
considered in the deconvolution procedure. 

3.5. Mean Instantaneous versus Effective Beam 

After applying the corrections described in Section 3.4, the 
resulting beam describes the telescope response to a point-like 
radiation source with approximately Rayleigh- Jeans (RJ) spec- 
trum, averaged over the focal plane of each array. We refer to 
this beam as the “mean instantaneous beam,” to emphasize that 
the combination of observations taken at many different times 
may lead to an effective window function that is different from 
the instantaneous one. 

For use in particular contexts, we compute an effective beam 
that includes corrections for alternative frequency spectra, and 
for pointing variation or other cumulative effects resulting from 
the combination of observations taken at many different times. 

We apply a simple first-order correction to obtain the effective 
beam for different radiation spectra. For radiation with a band 
effective frequency v, the beam is taken to be 

B\I) = B(£v rj/v), (12) 

where vrj is the effective frequency for radiation with an RJ 
spectrum. Effective frequencies for various spectral types are 
presented in Swetz et al. (201 1). 

In the absence of planetary sources, the absolute pointing 
registration of individual ACT observations is not known. Since 
each pixel of the ACT survey maps contains contributions from 
data acquired on many different nights throughout each season, 
the resulting season effective beam for these maps is less sharp 
than the instantaneous beam obtained from carefully aligned 
individual planet observations. 

Pointing repeatability can be estimated from the variation 
in apparent planet positions, relative to expectations, in each 
season and array. As summarized in Table 3, these indicate that 
the repeatability of telescope pointing is at the a p i ane t = 3" to 


8" (rms) level. Pointing variation in CMB observations may in 
principle be smaller than Opi ane t, since CMB observations are 
performed at the same azimuth angles each night, while planets 
observations span a wider range of azimuths. 

In addition to pointing repeatability, the season effective beam 
may also be diluted due to errors in the global pointing model, 
which is used to combine observations made at different azimuth 
angles. These global adjustments are estimated using point 
source positions measured in (non-cross-linked) maps made 
from only rising or only setting data. 

While pointing variance and global alignment may contribute 
to dilution of the season beam, we must also acknowledge that 
the beam may change slightly over the course of the season in 
a way that is not captured by the Saturn observations. (Changes 
in the beam over the course of the season, and the resulting 
impact on calibration, are assessed using Uranus observations in 
Section 4.) We thus seek an empirical measure of the difference 
between the season effective beam and the instantaneous beam 
obtained from Saturn observations. 

We parameterize the difference between the season effective 
beam Bf f and the instantaneous beam Bi as a Gaussian in l\ 

Bf = B t x e~^ +1 ^ 2 . (13) 

If the correction is interpreted as arising purely from Gaussian 
pointing error, then V is simply the pointing variance in 
square radians. While this motivates the form of the correction, 
our broader interpretation permits V to differ from the value 
expected based on pointing variance measured from planets, or 
to be less than zero if the season effective beam is somewhat 
sharper than the ensemble of Saturn observations indicates. 

We measure V for each season and array using bright point 
sources in the season CMB survey maps. Point sources with 
signal to noise ratios greater than 10 are identified using a 
matched filter. A beam model is fit to the map in the vicinity 
of each source, producing a value of V. The weighted mean 
and error of these individual fits, for each season and array, are 
used to form the effective beam for use with the survey maps. 
While this point source analysis gives us a secondary probe of 
the beam core, we note that the Saturn observations are still key 
for measuring the behavior in the wings of the beam. 

The values of the season correction parameter V are presented 
in Table 3. The measurement of V is somewhat susceptible to 
changes in the fitting parameters and thus the quoted errors 
have been inflated from the formal values by factors of two and 
six for the 148 GHz and 218 GHz arrays, respectively. While 
values of V for the 148 GHz array are very similar between 
seasons, values for 218 GHz are more variable. The analysis 
of Uranus observations in Section 4.1 shows that in 2010 the 
telescope beam was more sharply focused, relative to the mean 
beam estimated from Saturn, during most of the observing 
season (see Section 4.1); this results in a negative value for V 
at 218 GHz. 

4. PLANET BRIGHTNESS MEASUREMENTS 

In this section we use observations of Uranus and Saturn to 
study the variation in the system gain due to changes in water 
vapor level and telescope focus. 

4.1. Planet Amplitudes 

We characterize the apparent brightness of a planet in a map 
by examining the ratio of the Fourier transform of the map 
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to expectations based on the beam and planet shape. While 
ultimately equivalent to procedures that involve measuring 
source peak heights in filtered real-space maps, the Fourier space 
approach emphasizes the role of spatial filtering in providing 
precise planet amplitude measurements. It also allows us to 
quantify differences between the planet shape and expectations 
based on the beam measurements in a physically meaningful 
way. 

Once again working in the flat sky limit, we let 0 represent a 
vector in the plane of the sky centered on the planet. The Fourier 
conjugate variable is denoted as i and is written in subscript. 
The signal from a planet is modeled as a circular disk of uniform 
RJ temperature T p . For each planet observation, we use the total 
disk solid angle Cl p to compute angular radius 0 P = (Q p /7t) 1/2 . 
The Fourier transform of the planetary signal has components 

T t = T P Q P , (14) 

where Q p = £l p x 0 p )/(l 0 P )) is the Fourier transform of 
a disk of radius 0 P . For Uranus, Cl p ~ Q p at the ACT angular 
scales; for Saturn, the finite disk size cannot be neglected, 
but the effects due to its oblateness are negligible. Planetary 
radii are taken at the 1 bar atmospheric surface, as reported 
in Archinal et al. (2011). The Saturn equatorial radius is 
60268 zb 4 km and the polar radius is 54364 ±10 km; the 
Uranus equatorial radius is 25559 ± 4 km and the polar radius is 
24973 ± 20 km. 

Adapting Equation (8) to the flat sky, the ACT map M(0) of 
the planetary signal will have Fourier components 


M t = G 0 B l T p Q. P i . (15) 

For each planet map, we compute the amplitude M i of the 
planet at multiple bins in i. The map is centered, apodized, and 
Fourier transformed. To suppress atmospheric contamination, 
which contributes faint horizontal streaks to the map, we mask 
Fourier components having \t x \ < 1700. 18 The remaining 
complex Fourier amplitudes are averaged inside annuli of width 
A i ~ 600 and the imaginary part is discarded. We then compute 
the ratio of Mi to Bi£l p (which has been corrected to account 
for the apodization and masking applied to obtain Mi). 

The result is a curve of planet amplitude measurements, 


Oil = 


Mi 

B t Cl p ’ 


(16) 


for 2000 < i < 14, 000. In the absence of any gain or beam 
variation, oti should be equal to GoT p , the apparent brightness 
of the planet in detector power units. 

The mean and variance of ai curves for selected seasons 
and arrays are shown in Figure 4. The ai curves are not in all 
cases consistent with a constant value. This may be attributed to 
variations in the focus from night to night. To parameterize these 
deviations and to make possible an assessment of the impact they 
might have on the system calibration, we fit parameters a and 
P for each ai curve according to the following model: 

ai=all + P . (17) 

V 5000/ 7 


18 No such masking or filtering is performed on the Saturn maps prior to their 
use in the determination of the telescope beam; the masking is more important 
for Uranus maps, because the signal is weaker. 


148 GHz - 2008 


218 GHz - 2010 



4000 


8000 

t 


12000 




Figure 4. Variation of planet amplitude measurements with angular scale. Each 
panel shows mean fractional deviation of ai relative to a (solid line) and the 
standard deviation of all individual at fa curves (gray band) for a particular 
season, array, and planet. See Section 4. 1 for definitions of ai and a . We 
show results from 148 GHz 2008 (left panels), where there is a high degree of 
consistency between the planet observations and the season mean beam. We also 
show results for 218 GHz 2010 (right panels), where there is somewhat more 
variance in the curves (which is typical of the 218 GHz array) and where the 
Uranus observations deviate significantly from the mean beam (which is seen 
in all arrays in 2010). Differences between Uranus and Saturn are indicative 
of 2% changes in the beam focus during the season; such changes have been 
accounted for in the season effective beam by fitting a correction parameter to 
maps of point sources in the season maps as described in Section 3.5. 


The mean angular scale i isjuken as the center of the fitting 
range. The parameters a and P represent the mean amplitude of 
the planet and the deviation Jfom expected focus of the beam, 
respectively. For example, P = 0.01 would indicate that the 
planetary signal is stronger by 1% at i ~ 10, 000 relative to 
i ~ 5000, given the mean season beam, and thus the telescope 
beam at the time of observation was slightly more sharply 
focused. We would expect observations associated with negative 
ft to also have a smaller a than expected, since a less sharply 
focused beam will have poorer overall efficiency. 

The extraction of amplitudes a and focus parameters p allows 
us to quantify deviations of focus over the season, and to estimate 
the effects of such deviations on the overall system calibration. 
This is accomplished in the next section, where a model for 
system calibration is obtained from the brightness and focus 
measurements of Uranus. 

4.2. Calibration Parameters from Uranus Observations 

The measurements of a and ft from Uranus observations, 
along with measurements of the precipitable water vapor (PWV) 
level for most of these observations, permit us to characterize 
the impact of PWV level and focus variation on the telescope 
calibration. This analysis does not require that the brightness 
temperature of Uranus be precisely known, but assumes that the 
brightness is roughly constant over the observing season. (We do 
not perform the same analysis on Saturn, because the effective 
brightness varies with the ring opening angle; see Section 5.2.) 

The Atmospheric Transmission at Microwaves (ATM) model 
(Pardo et al. 2001) provides estimates of atmospheric emission 
and absorption at millimeter and sub-millimeter wavelengths. 
The ACT frequency bands avoid strong molecular resonances 
and are thus primarily susceptible to continuum emission and 
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absorption due to water vapor. For the purposes of calibration 
we are interested in the atmospheric transmission in a given 
frequency band. The transmission at the ACT observation 
altitude Go is written as 

T = exp(— (r w u; + Td)/sin0o), (18) 

where w is the PWV column density at zenith and z w w and r d 
are the “wet” and “dry” optical depths at zenith. The parameters 
t w and id vary negligibly over the range of temperatures and 
pressures experienced at the telescope site. 

Because planet observations discussed in this work are ob- 
tained at the same fiducial altitude used for CMB observations, 
the effects of Td are common to both observations and cannot be 
measured from the planet observations. We define the overall 
system gain Go such that it applies at the typical PWV level of 
w 0 = 0.44 mm. The w -dependent gain factor is then 

Gq(w) = G 0 exp (— (t w + z x )(w - wo)/ sin 0 ) . (19) 

Here r w is as provided by the ATM model, while t x is a fit 
parameter that could indicate a deviation from the ATM model 
but is more likely due to a bias in the detector calibration 
procedure that is sensitive to sky loading. Estimates of r w for the 
ACT bands were computed using the ATM model by Marriage 
(2006); these values are provided for reference in Table 4. 

To assess the impact of small variations in telescope focus on 
planet calibration, we include jt dependence on the difference 
between the focus parameter /3 (defined in Equation (17)) and 
the season mean focus parameter fi. Because a planet was 
observed approximately once per night, we take /3 to be the 
average of measured /3 for all planet observations for the array 
and season under consideration. We adopt an exponential form 
out of convenience; calibration variation is small so a linear 
correction would be equivalent. 

The full model for the Uranus brightness data a, is then 
given by 

a = T u G 0 e m ^~^ exp (— (r w + z x )(w — wo)/ sin0) , (20) 

where the fit parameters T u Go, m, and r x represent overall 
calibration, sensitivity to variations in telescope focus, and sen- 
sitivity to PWV beyond predictions based on ATM, respectively. 

This model is fitted to the Uranus brightness amplitude data 
for each season and frequency array. The parameter values are 
presented in Table 4. The residuals have scatter of approximately 
2%, 2%-6%, and 6% in the 148 GHz, 218 GHz and 277 GHz 
arrays respectively; this scatter is not explained by uncertainties 
in planet amplitude measurement. Values obtained for r x are in 
good agreement between seasons for the 148 GHz and 218 GHz 
arrays. They are inconsistent with zero (at the ~2cr to 4cr level), 
which we attribute to systematic detector calibration error. For 
the 277 GHz array there is some slight disagreement between the 
2009 and 2010 results, which we attribute to the low number of 
observations in 2009. The values of t x are used in combination 
with r w to correct the time-ordered data for opacity prior to 
creation of full-season survey maps, as described in Diinner 
et al. (2013). 

While for some seasons and arrays we obtain better fits to the 
data by including the parameter m, the impact of this additional 
parameter on inferred overall system calibration is small. To 
show this, we refit the model of Equation (20) with parameter 
m fixed to 0. The resulting calibration factor T u G^^ may be 
compared to the case where m was free to vary. The differences 


Table 4 

Properties of Planetary Data and Calibration Model 
Fit Parameters by Season and Array 




Array 


Quantity 

148 GHz 

218 GHz 

277 GHz 

r w (mm -1 ) 

0.019 

0.044 

0.075 

r x (mm -1 ) 




2007 

0.016 ± 0.010 



2008 

0.014 ± 0.004 

0.040 ± 0.008 

0.033 ± 0.063 

2009 

0.009 ± 0.005 

0.021 ± 0.017 

-0.008 ± 0.041 

2010 

0.012 ± 0.003 

0.045 ± 0.008 

0.083 ± 0.016 

Error inT^Go (%) 



2007 

0.9 



2008 

0.2 

0.4 

1.3 

2009 

0.3 

u 

2.2 

2010 

0.2 

0.4 

0.8 

m 




2007 

0.7 ± 0.6 



2008 

0.5 ± 0.3 

-0.5 ± 0.6 

1.2 ± 1.3 

2009 

0.5 ± 0.3 

3.4 ± 0.5 

1.2 ±0.7 

2010 

0.8 ±0.1 

0.5 ± 0.2 

1.5 ±0.4 

/3 mean (rms) 



2007 

-0.017 (0.010) 



2008 

0.002 (0.006) 

0.011 (0.010) 

0.011 (0.015) 

2009 

0.004 (0.011) 

0.008 (0.016) 

0.007 (0.020) 

2010 

0.015 (0.015) 

0.021 (0.019) 

0.017 (0.019) 

(G^= 0) /G 0 - 

- 1) (%) 



2007 

-0.6 



2008 

0.1 

-0.3 

0.5 

2009 

0.2 

1.4 

0.4 

2010 

0.3 

0.2 

0.4 


Notes. The opacity parameters r w from the ATM model are 
provided for reference. Fit parameters r x and m parameterize 
sensitivity of Uranus apparent brightness to PWV and focus 
parameter, respectively. The mean and standard deviation of 
the focus parameter /3 for all planet observations in each 
season is provided for reference. The overall calibration shift 
G ( ™- 0) /G 0 — 1 is the systematic calibration change if m is fixed 
to zero instead of being free to vary. 

are presented in Table 4; the change in overall calibration is 
0.6% or smaller, and always less than the la error on T u Go, 
with the exception of the 2009 data for the 218 GHz array. This 
exceptional case also has a large (9%) uncertainty in its absolute 
WMAP - based calibration. 

The fit values of T u Go may be used to provide an absolute 
calibration of each array (should T u be known) or to provide 
a measurement of the Uranus disk temperature based on an 
independent calibration of Go. The latter possibility is discussed 
in the next section, leading to the measurements of Uranus and 
Saturn brightness temperatures presented in Section 5. 

4.3. Relation to WMAP -based Calibration 

The absolute calibration of the ACT full- season CMB survey 
maps is obtained through cross-correlation in £-space (Hajian 
et al. 2011) to the 94 GHz maps from WMAP1 (Jarosik et al. 
2011). As described by Das et al. (2013), the WMAP maps are 
used to calibrate the ACT 2008 Southern and 2010 Equatorial 
field maps at 148 GHz directly; the latter maps are then used to 
calibrate the 148 GHz map for 2009 and the 218 GHz maps for 
all seasons. 
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While the data entering the survey maps are corrected for 
atmospheric opacity, they have not been corrected for any kind of 
time variation in the telescope beam. Instead, we have obtained 
an effective season beam suitable for use with the survey maps 
by comparing our mean instantaneous beam to point sources in 
the survey maps, as described in Section 3.5. This correction has 
negligible effect at low i where the calibration to WMAP takes 
place. Furthermore, the WMAP calibration should be associated 
with the telescope gain at the season- average focus parameter. 
Thus the calibration factor obtained from the WMAP calibration 
is compatible with Go as defined by Equation (20). 

5. BRIGHTNESS TEMPERATURES OF 
URANUS AND SATURN 

Using the WMAP - based calibration of the ACT survey maps 
at 148 GHz and 218 GHz, we convert the Saturn and Uranus 
amplitude measurements to calibrated brightness temperatures. 
A geometrically detailed, two-parameter model for Saturn disk 
and ring temperatures is fit to the Saturn brightness data. Because 
CMB survey maps for the 277 GHz array have not yet been 
calibrated to WMAP , we do not consider that array in this 
analysis. 

In what follows, we refer to both the RJ and brightness 
temperatures of the planet. 19 The RJ temperature T RJ is related to 
the brightness temperature by equating the specific intensity 
of an RJ spectrum and a blackbody spectrum at the band effective 
frequency v e ff! 

B v JT b ) = 2v 2 efl k B T R1 /c 2 , (21) 

where B V (T ) is the blackbody spectral radiance. The spectra of 
Saturn and Uranus are each sufficiently close to the RJ limit 
that we take the source effective frequencies to be the RJ band 
centers of 149.0 GHz and 219.1 GHz (Swetz et al. 2011). 

5.1. Uranus 

Based on our fitted values of T u Gq (Section 4.2) and 
the absolute calibration of Go using WMAP data, we obtain 
measurements of T u . Because the calibration factor is obtained 
from the primary anisotropies of the CMB, we must account for 
the different frequency spectrum of the planets relative to the 
CMB blackbody. The effective spectral index of the flux of the 
CMB in the ACT bands is 1.0 (0.0) in the 148 GHz (218 GHz) 
bands, while the spectral index of Uranus is approximately 
1.65 in both bands (based on, e.g., Griffin & Orton 1993). The 
3.5 GHz uncertainty in the ACT band centers and the effective 
frequencies provided in Swetz et al. (201 1 ) result in an additional 
1.6% (2.7%) error in the Uranus brightness determinations at 
148 GHz (218 GHz). 

The RJ temperatures obtained for each season and array 
are presented in Table 5. The value quoted for the 148 GHz 
(218 GHz) array is valid at effective frequency of 149.0 
(219.1) GHz. All RJ temperatures have been inflated by 0.5 
(0.3) K to compensate for the background provided by the CMB 
in our bands. Brightness temperatures may be computed from 
the RJ values by adding 3.5 (5.1) K. 

The results are consistent between seasons (even after re- 
moving the correlated error due to the frequency spectrum cor- 
rection), with uncertainty dominated by the overall instrument 

19 Our planet amplitude measurements include a deconvolution of the 
telescope beam and thus we do not work with the planet’s “antenna 
temperature” directly. 


Table 5 

RJ Temperature Measurements of Uranus and Saturn 



7* J (JQ 

148 GHz 218 GHz 

Uranus (Section 5.1) 

2008 

101.5 ±2.0 

93.2 ± 2.3 

2009 

103.8 ± 2.2 

91.5 ±8.6 

2010 

105.7 ± 2.1 

97.6 ± 2.6 

Combined 

103.2 ± 2.2 

95.0 ±3.1 

Saturn, Wll model (Section 5.2) 

Tlisk 

133.8 ±3.2 

132.2 ±4.7 

Ting 

17.7 ± 2.2 

12.3 ± 4.0 

[x 2 / d.o.f.] 

[7.6/5] 

[6.5/5] 

Saturn, single T model (Section 5.2) 

Tlisk 

131.2 ±3.1 

127.5 ± 4.4 

[X 2 /d.o.f.] 

[12.0/6] 

[11.2/6] 


Notes. Values for Uranus (Section 5.1) are shown for each season, as 
well as a weighted mean that takes full account of calibration covariance 
between seasons. Saturn results (Section 5.2) are provided for both a 
two-component disk+ring model and a disk-only model. All temper- 
atures are RJ and have been corrected for the CMB to indicate the 
brightness that would be seen relative to an empty (T = 0) sky. To 
obtain thermodynamic brightness temperatures at 148 GHz (218 GHz), 
add 3.5 K (5.1 K) to the RJ values. 


calibration uncertainty rather than statistical uncertainty from 
planet brightness measurements. 

From the RJ temperature measurements in each season, we 
compute a weighted mean disk temperature for Uranus. The 
WMAP - based calibration error is correlated between the 2009 
and 2010 seasons, because the 2009 map is calibrated to the 2010 
ACT map, and the ACT band center uncertainty is correlated 
across all seasons. We account for this in the mean and error 
presented. At 148 GHz, we obtain a RJ (brightness) temperature 
of 103.2 K (106.7 K), with 2.1% error. At 218 GHz, we obtain 
temperature 95.0 K (100.1 K), with 3.2% error. 

We compare our Uranus temperature measurements to the 
widely used model of Griffin & Orton (1993, hereafter G&O). 
The G&O model provides brightness temperatures for Uranus 
and Neptune in the ~100 to 1000 GHz range based on precise 
measurements of brightness ratios of each planet to Mars. The 
absolute calibration is obtained by interpolating, in the logarithm 
of the frequency, between the Wright (1976) model for Mars 
temperature at 3.5 /zm (857 GHz), and the (Ulich 1981) model 
for Mars temperature at 90 GHz. They adopt a 5% systematic 
error for this Ulich-Wright hybrid model. 

The G&O model predicts a brightness temperature of 112 K 
at 148 GHz, approximately 5% larger than our measured value. 
This suggests that the Ulich-Wright model for Mars brightness 
is 5% high at 148 GHz. Weiland et al. (2011) have observed 
a similar 5% discrepancy between Mars temperatures and the 
Ulich model at 94 GHz. At 218 GHz, the G&O model gives a 
brightness temperature of 99 K, which is 1 % smaller than our 
measured value; this difference is smaller than our calibration 
uncertainty. 

Long term studies of Uranus show a variation in mean disk 
brightness at 8.6 GHz (Klein & Hofstadter 2006) and 90 GHz 
(Kramer et al. 2008) on decade time scales. If this variation 
is attributed to differences in the Uranus surface brightness 
at different planet latitudes, it implies that the mean disk 
temperature is roughly 10% larger when the South pole, rather 
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than the equator, is directed to toward the observer. 20 The ACT 
measurements over the 2008-2010 period span Uranus sub- 
Earth latitudes from 1° to 13°. In comparison, the observations 
presented by G&O were made in 1990-1992 and span sub-Earth 
latitudes from —70° to —60°. It is then possible that, rather than 
a difference in calibration, Uranus was brighter in these bands at 
the time the G&O data were taken. However, if one accepts the 
5% downward recalibration of the Ulich Mars model at 94 GHz, 
then one should accept a downward recalibration (of 4%-5%) 
of the G&O data around 148 GHz. This brings the ACT and 
G&O measurements into agreement, with no evidence for any 
variation in Uranus brightness temperature due to sub-Earth 
latitude. 

Our results at 148 GHz are consistent with the analysis of 
Sayers et al. (2012), who use WMAP measurements of Uranus 
and Neptune at 94 GHz (Weiland et al. 2011) to recalibrate the 
G&O model at 143 GHz. It is important to note that while our 
measurements rely on a detailed understanding of the telescope 
beam, our absolute calibration standard is the CMB, for which 
the brightness as a function of frequency is known. Our errors 
thus do not include any significant systematic contribution due 
to interpolation of unknown source spectra into our frequency 
bands. 

In Figure 5 we show our brightness temperature measure- 
ments along with WMAP measurements below 100 GHz from 
Weiland et al. (2011) and with data used by G&O to fit their 
model in the 90-1000 GHz range. 

As discussed above, systematic error in the G&O model 
is dominated by uncertainty in the Mars spectrum below 
1000 GHz. The ACT and WMAP data provide new, absolutely 
calibrated temperature measurements between 20 and 220 GHz. 
We thus fit a simple empirical model for the temperature using 
the ACT, WMAP , and the higher frequency G&O data. To avoid 
the strong spectral feature at 30 GHz, we include only the two 
highest frequency bands of WMAP. We use the three data points 
in G&O above 600 GHz. We model the Uranus brightness 
temperature as a function of 0 = (v/100 GHz) by 

'j'U 

— =a o + ailog lo </> + a2logi O 0, (22) 

1 JV 

and obtain best- fit coefficients (« a o, a\, < 22 ) = (121, —78.2, 18.2). 

The ACT, WMAP , and G&O brightness temperature mea- 
surements are shown in Figure 5, along with the G&O model 
and our empirical model. In this comparison, we include the 
ESA2 and ESA4 models, which were brought to our attention 
after the initial submission of this work. 21 

5.2. Saturn 

Measurements of the Saturn disk temperature are complicated 
by the presence of the rings, whose dust both obscures the main 
disk and radiates at a lower temperature. The ACT observations 
of Saturn span ring opening angles from —2° to 12° and provide 
an opportunity to explore the separate contribution of these two 
components to the total effective brightness of the planet. 

We follow closely the approach of Weiland et al. (2011, 
hereafter Wll), who apply a two-parameter model in which 


20 The center of the projected planet disk, as observed from Earth, is called the 
sub-Earth point. The Uranian latitude of this point is called the sub-Earth 
latitude. 

21 Models retrieved on 2013 June 18 from 

ftp://ftp.sciops.esa.int/pub/hsc-calibration/PlanetaryModels/ESA2/ and 
ftp://ftp.sciops.esa.int/pub/hsc-calibration/PlanetaryModels/ESA4/. 
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Figure 5. Measurements of the Uranus brightness temperature from this work, 
Weiland et al. (2011, WMAP), and Griffin & Orton (1993, G&O). The G&O 
points include data from that work as well as Ulich (1981) and Orton et al. 
(1986). The dotted line shows the model of G&O. G&O data are calibrated to 
a Ulich-Wright hybrid model for Mars brightness that interpolates between 90 
and 857 GHz and carries an estimated 5% systematic error. The solid line is 
our best- fit empirical model, using only the two highest frequency points from 
WMAP, the two ACT points, and the G&O points above 600 GHz. Fractional 
difference of data and models relative to our best-fit empirical model is shown 
in bottom panel. 

(A color version of this figure is available in the online journal.) 


the opacity and geometry of the rings are fixed, while the disk 
temperature and an effective temperature for the rings are free to 
vary. The ring dimensions and normal opacities are taken from 
Dunn et al. (2002). For each of the 148 GHz and 218 GHz arrays, 
we fit Tdisk and Turing ? which describe the mean disk temperature 
and mean effective brightness temperature of the rings (see Wll 
for the detailed definition), respectively. 

The measurements of a and for each Saturn observation 
are converted to RJ temperatures, including corrections for PWV 
level and focus parameter, and with the CMB -based calibration 
applied. Because the scatter in the temperatures is not fully 
accounted for by the errors in the amplitude measurements and 
correction parameters, we combine observations made in single 
15-30 day periods, taking the uncertainty in each combined 
measurement to be the error in the mean of the contributing data. 
As a result, a small number of temporally isolated observations 
are discarded. The binned data points, and the best-fit model, 
are shown in Figure 6. 

The fit of the Wll model includes a full accounting of the 
non-trivial correlations in the calibration uncertainty of the data 
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Figure 6. Binned data and resulting best-fit two-component model of effective 
Saturn brightness vs. ring opening angle ( B ) relative to observer. The model is 
symmetric about B = 0 by construction. As the absolute value of B increases, 
the increased radiation from the rings and decreased flux from the obscured disk 
leads to a local minimum in the total effective brightness at B ^9° (13°) at 
149 GHz (219 GHz). The model is fitted independently for each frequency band. 
Error bars correspond to the error in the mean of observations contributing to 
each point, but fits include additional error due to calibration, which is covariant 
within each season and in some cases between seasons. Approximate mean 
dates of observation are, from left to right: 2008 November, 2008 December, 
2010 June, 2010 April, 2009 December, 2010 December, 2010 Decemeber. 
The sign convention is such that B < 0 corresponds to negative values of the 
sub-Earth latitude. 
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Figure 7. Measurements of the Saturn disk brightness temperature from ACT, 
Weiland et al. (2011, WMAP ), and Goldin et al. (1997, G97). Points from G97 
have been recalibrated as described in the text. Frequency error bars on ACT 
and G97 points indicate bandwidth. Temperature error bars on G97 points do 
not include calibration uncertainty. The green dashed line is a spectral model 
from Weisstein & Serabyn (1996, WS96), which shows proximity of absorption 
features to G97 points. 

(A color version of this figure is available in the online journal.) 
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points. The resulting models are good fits, with x 2 per degree of 
freedom of 7.6/5 (6.5/5) in the 148 GHz (218 GHz) array. The 
model fit parameters and uncertainties are presented in Table 5. 

For comparison, we also fit a single temperature model to the 
binned Saturn data, finding poorer fits and mean temperatures 
that are 1 a less than the disk temperatures obtained in the full 


model. These temperatures and x 2 statistics are also presented 
in Table 5. 

In Figure 7 we show the ACT disk temperature measurements 
along with high precision results from Wll below 100 GHz 
and from Goldin et al. (1997, hereafter G97) between 170 
and 700 GHz. The G97 measurements have statistical errors 
at the 2 % level, but are calibrated with reference to the 
same Mars model used by G&O. We have applied a rough 
recalibration factor to the G97 measurements using the WMAP 
result (Weiland et al. 2011) that the Ulich (1981) model for 
Mars temperature is 5% high at 94 GHz. Our recalibration 
factor varies linearly with the log of frequency from 0.95 at 
90 GHz to 1 at 857 GHz. The brightness temperature at 170 GHz 
presented by G97 lies somewhat above the ACT measurements 
at 149 and 219 GHz. However, with our recalibration factor 
applied, the G97 and ACT data are consistent with a constant 
brightness temperature over this frequency interval. It is difficult 
to comment further on the continuum spectrum of Saturn based 
on the ACT and G97 data, due to the proximity of the high 
frequency G97 bands to PH 3 resonances (e.g., Weisstein & 
Serabyn 1996). 


6 . CONCLUSION 

We have described the reduction of ACT planet observations 
for the purposes of measuring and monitoring the telescope 
beam and gain parameters. The core of the beam has been 
well characterized by Saturn observations, and studies of point 
sources in the map have provided an empirical correction 
to account for pointing variance, alignment error, and slight 
changes in beam focus. Observations of Uranus have been 
used to measure and account for changes in system gain due 
to atmospheric water vapor and focus variation. 

Based on hundreds of brightness measurements of Uranus, 
we have obtained precise measurements of the disk bright- 
ness temperature at 148 GHz and 218 GHz. In contrast to 
previous observations, our absolute calibration standard is the 
CMB blackbody. We have also obtained precise measurements 
of the Saturn disk temperature, in the context of a simple 
two-component model. 
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